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ABSTRACT 



We study the frequency dependence of the frequency shifts of the low-degree p modes measured in the F5V star HD 49933, by 
analyzing the second run of observations collected by the CoRoT satellite. The 137-day light curve is divided into two subseries 
corresponding to periods of low and high stellar activity. The activity-frequency relationship is obtained independently from the 
analysis of the mode frequencies extracted by both a local and a global peak-fitting analyses, and from a cross-correlation technique 
in the frequency range between 1450 //Hz and 2500 pHz. The three methods return consistent results. We show that the frequency 
shifts measured in HD 49933 present a frequency dependence with a clear increase with frequency, reaching a maximal shift of about 
2 /jHz around 2100 //Hz. Similar variations are obtained between the / = and 1=1 modes. At higher frequencies, the frequency 
shifts show indications of a downturn followed by an upturn, consistent between the / = and 1 modes. We show that the frequency 
variation of the p-mode frequency shifts of the solar-like oscillating star HD 49933 has a comparable shape to the one observed in the 
Sun, which is understood to arise from changes in the outer layers due to its magnetic activity. 

Key words. Stars: oscillations - Stars: solar-type - Stars: activity - Methods: data analysis 



1. Introduction 

The convective motions in the external layers of solar- 
like oscillating stars excite acousti c sound waves, which are 
trapped in the stellar interio rs (e.g., iGoldreich & Keelevl 1 19771 
iGoldreich & Kuman, 119881) . The precise frequencies of these 
pressure (p) driven waves depend on the properties of the 
medium in which they propagate. Thus, helio- and asteroseimol- 
ogy are able to infer the properties of the Sun and stellar 
interiors by studying and characterizing t hese p modes (e.g., 
iGough et all 1 19961: iMetcalfe et all l2010bl) . In particular, since 
the very first helioseismic observations, different authors no- 
ticed that near the maximum of the solar magnetic activity cycle, 
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19921). The frequency shifts, which 



follow an inverse mode-mass scaling dLibbrecht & Woodardl 
1990), are explained to arise from cha nges in the outer lay- 
ers of the Sun along the 1 1 -year cycl e (IGoldreich et al., 1991; 
iDziembowski & GoodeL I2004L l2005h . The high-frequency p 
modes, which have higher upper turning points, are then more 
sensitive to the perturbations induced by the magnetic field in 
the outer part of the Sun than the low-frequency p modes, that 
are less sensitive to the outer layers. Thus, the high-frequency p 
modes are observed to have large variations with activity, while 
the low-freque ncy p modes present little or no variation with 



As longer, continuous, high-quality helioseismic measure- 
ments were collected, the other p-mode oscillation parame- 
ters, such as the amplitude, the linewidth, and the asymmetry, 
were also observed t o vary with the sol ar activity cycle (e.g , 
[ Chaplin et ail l2000t iKomm et all 120021: ISalabert et all 12003 ; 
Salabert & Jimenez-Reves, 2006; Jimenez-Reves et al., 2007). 
While the frequency shifts were proven to be closely correlated 
with the surface activity proxies over the las t three solar cycles 
(21, 22, and 23) (e.g.. IChaplin et all l2007bl) . clear differences 
were recently shown between the global p-mode frequencies and 
the surface activity of t he Sun during the unu s ually extended 
minimum of cycle 23 dBroomhall et all 120091: Salabert et al., 
2009; iTripathv et all l2010h . Furthermore, a faster cyclic fre- 
quency variations, with a period of about 2 year s , was observed 
to co exist in the Sun (Benevolenskayal 119951 : iFletcher et all 
1201 Oh . which is believed to have an origin located deeper inside 
the convective zone. 

Magnetic cycle s in othe r stars were already reported by sev- 
IWilsonl 11978b iBaliunas & Vaughanl Il985t 
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activity (e.g., Libbrecht & Woodard, 1990; Ronanetal., 
Chapli n et all 119981: iGellv et all 120021: ISalabert et all 
Howe et alll2008l) . 
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eral authors 

IBaliunas et all 119951; lHall. Lockwood & Skiffl. 12007b . and the 
observed cycles are covering a range between 2.5 and 25 years. 
Moreover, several authors suggested that the periods of the activ- 
ity cycles increase proportional to the stellar rotational periods 
along two distinct paths in main-sequence star s : the active and 
the inactive stars (Saar & Brandenburg, 119991 iBohm-Vitensel 
120071) . Recently. IMetcalfe et all d2010at) uncovered a 1.6-year 
magnetic activity cycle in the exoplanet host star HD 17051 by 
measuring the Ca n H and K emission lines. However, techniques 
based on spectroscopic measurements are only good proxies of 
surface magnetic fields (LeightorJ Il959l) . Spectropolarimetric 
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observations allow direct measurements of stellar magnetic 
fields, that has revealed, for instance, reversals of magnetic field 
polarity during stellar cycles for F and G stars (see lDonati et al.L 
2008: lFares etafl 120091: IPetit et al.Ll2009h . 

The recent discovery of variations with magnetic activity in 
the p-mode oscillation frequenci es of the F5V star H D 49933 
using asteroseismic techniques (Gar cia et all [2010) opens a 
new era in the study of the physical phenomena involved in 
the dynamo processes. Indeed, the seismic observables pro- 
vide invaluable information to pierce inside the stars and de- 
termine key parameters for the study of stellar activity, such 
as the depth of the convection zone, the characteristic evolu- 
tion time of the granulation, the differential rotation, or the 
sound-speed as a function of the star's radius. All these new 
observables will impose new constraints to the theory (e.g. 
IChaplin et all l2007at iMetcalfe et all 120071) . Long, high-quality 
photometric observations being currently collected by the space- 
based Convection, Rotation, a nd planetary Tra nsits (CoRoT, 
iBaglin et alll2006l) and Kepler (iKoch et allboiOl) missions will 
certainly contribute to these studies (e .g. Barb an et al. 1 120091 
Garcia et all 120091; iDeiieuvels et alll201ollMafhur et all 1201 Oat 
Chaplin et al.Ll2010h. In this p aper, we extend the original anal- 



ysis from lGarcfa et al. (2010) by studying the frequency depen- 
dence of the activity-frequency relationship. 

2. Data analysis 

The light curve of HD 49933 was col lected during two observa- 
tional runs from the C oRoT satellite ( Appourch aux et al.l [2008; 
Benomar et al ., 2009). The first run was performed in 2007 over 
a period of 60 days, while in 2008 a second run of observations 
was obtained over 137 days. In the present paper, we analyzed 
the data from the second run only as it is the longest run and 
that it corresponds to an important variation in the magnetic ac- 
tivity of HP 49933 and its p-mode frequencies and amplitudes 
(IGarcfa et al.L l2010). These continuous 137 days of observations 
were divided into two subseri es depending on the activity level 
of HD 49933 as observed by IGarcfa et all (120101) . The first 90 
days of the original time series correspond to a period of low ac- 
tivity, and the following 47 days correspond to a period of higher 
activity. Thus, these two subseries optimize the temporal cover- 
age of the observed magnetic activity as well as the frequency 
resolution. The frequency dependence of the frequency shifts of 
HD 49933 were obtained by three independent methods: 1) by 
means of the mode frequencies extracted by a local (method #1) 
and global (method #2) peak-fitting procedures, and 2) by using 
a cross-correlation analysis (method #3). Details of the methods 
are given below. A total of twelve overtones over the frequency 
range from 1450 /vHz up to 2500 /iHz of the oscillation power 
spectrum were analyzed. 



2.1. Mode peak-fitting methods 

The power spectrum of the two subseries was obtained in or- 
der to extract estimates o f the p-mode p a ramet ers. The mode 
identification provided by iBenomar et al.l ({2009) was adopted. 
The stellar background wa s modeled using six free parameters 
as in Mathur e t ak| (|2010bl) and the individual frequencies from 
IBenomar et al.l (I2009I) were taken as initial guesses. Both local 
and global peak-fitting techniques were used. The local peak- 
fitting procedure was performed by fitting sequences of succes- 
sive series of I = 0, 1, and 2 modes using a maximum-likelihood 
estimator (method #1). The individual / = 0, 1, and 2 modes 
were modeled using a single Lorentzian profile for each of the 



angular degrees I. Therefore, neither rotational splitting nor in- 
clination angle were included in the fitting model. The amplitude 
ratios between the / = 0, 1, and 2 modes were fixed to 1, 1.5, and 
0.5 respectively, and only one linewidth was fitted per radial or- 
der. Although the I = 2 modes are visible in the power spectrum, 
the accuracy of their frequency measurements is rather small 
co mpared to the / _ and 1 modes (see for instance, Table 1 
of Benomar et al. , 2009). We decided then to use only the I — 
and 1 modes to calculate the frequency shifts from the individ- 
ual mode frequencies extracted with the method #1. Because the 
very short mode lifetimes in HD 49933, the I = 2 modes are 
blended with the neighboring I = modes, which makes the 
determination of the / = and 2 mode frequencies strongly neg- 
atively correlated, especially at high frequency where the modes 
are broader. Due to this strong blending, we decided to determine 
the centroid of the even I modes by fitting each pair / = and 2 
with a single Lorentzian profile. We then performed a global fit- 
ting of the spectrum including the stellar background and about 
20 Lorentzian profiles (method #2). The frequencies obtained for 
the even pairs of modes correspond to a weighted average of the 
I = and I = 2 frequencies, the weights being the relative visi- 
bilities of the modes. If we assume: 1) that the relative amplitude 
between the 1 = and 2 modes does not change either with fre- 
quency or with the activity level, and 2) that the frequency shifts 
depend mainly on the frequency, i.e. that close pairs of I = 
and 2 modes have similar activity-induced frequency variations, 
then the shift measured for the even centroids is equal to the shift 
of the / = modes. Note also that the fitted frequencies for the 
odd pairs of modes correspond to the I = 1 frequencies, as the 
I = 3 modes have neglecting contributions. The mode parame- 
ters thus extracted from method s #1 and #2 were chec ked to be 
consistent within the errors with IBenomar et ai1 (f2009). For both 
methods #1 and #2, the frequency shifts 5vi were obtained by 
subtracting the individual frequencies measured during the 47- 
day period of high magnetic activity and the 90-day period of 
low magnetic activity (similarly to what is performed with the 
Sun, see e.g., ISalabert et al.L 120041) . Additionally, we analyzed 
the temporal variations of the mode amplitudes and linewidths, 
but the results are not reliable due to a poor determination of 
these parameters at the level of precision of the data. 

2.2. Cross-correlation method 



The cross-correlation method was first used bv lPalle et al.l (fl989) 
in order to measure the solar frequency shift from one single in- 
strument. In our case, the power spectrum of the 47 days at high 
activity was correlated with the power spectrum of the 90 days 
at low activity (method #3). As the two datasets do not have the 
same number of days, the shortest subseries was zero-padded in 
order to have the same number of bins in both power spectra. 
The background noise was modeled as in Sect. 12.11 and sub- 
tracted from the power spectra. To obtain the frequency shift, a 
Gaussian profile was fitted to the cross-correlation function to 
estimate its maximum, the Gaussian profile being centered at 
the position given by the skewness of the cross-correlation func- 
tion to avoid any bias due to asymmetry. Thus, the centroid of 
the fitted Gaussian returned the amount of shift in frequency, 
the associated uncertainty being given by the lcr error bar of 
the fit. While the local and global peak-fitting methods allow 
us to study the frequency shifts of the individual modes sepa- 
rately (Sect. 12. U . the cross-correlation method returns a global 
value of the frequency shift for all the visible modes in the power 
spectrum. In order to obtain estimates of the frequency shifts for 
the even and odd pairs of modes, two set of spectra were ob- 
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Fig. 1. Smoothed power spectra of the 90-day period of low stellar activity (black) and of the 47-day period of high activity (red) of 
the analyzed observations of HD 49933 represented over the frequency range between 1440 and 2170 //Hz. 



tained from each spectra. These new spectra were built by ex- 
tracting the points containing only the frequency windows of the 
corresponding odd or even pairs, and creating thus a new array 
of points. Then, the cross-correlation functions were calculated 
in the same manner as previously and the associated frequency 
shifts of the even and odd pairs of modes returned. 

3. Results 

The power spectra of HD 49933 corresponding to the 90-day 
period of low stellar activity and the 47-day period of higher 
activity are shown on Fig. Q] between 1440 and 2170 //Hz. For 
illustrative purpose, the spectra were heavily smoothed over 20 
//Hz, which is compara ble to twice the mod e linewidth at high 
frequency (see Fig. 4 in lBenomar et all 2009). Both power spec- 
tra overlap in frequency in the low-frequency range below 1780 
//Hz, while a positive shift becomes visible as the mode fre- 
quency increases during the period of higher stellar activity. We 
note the presence of a prominent peak at about 1690 //Hz in the 
47 -day power spectrum of high activity. We checked that it does 
not correspond to a harmonic of the spacecraft orbit, therefore 
it might have an instrumental origin. Figure |2] shows the fre- 
quency shifts of HD 49933 as a function of frequency measured 
between the periods of high and low activity using the results 
from the local peak-fitting analysis (method #1). Similar results 
are obtained with the global peak-fitting analysis (method #2). 
A frequency dependence of the / = and 1 p-mode frequency 
shifts can be observed. The mean frequency shifts between con- 
secutive I = and 1, using the formal errors as weights, are 
also represented. The frequency shift increases with frequency 
up to about 2100 //Hz to reach a maximum of ~ 2 /iHz, followed 
by a downturn with a reversal of the sign down to about —1.5 
//Hz at ~ 2300 //Hz for the mean frequency shift. Moreover, 
this downturn might be followed by an upturn for frequencies 
above 2300 //Hz. The individual I = and I = 1 modes have, 
within the errors, similar variations with frequency. But the / = 1 
frequency shifts are overall better constrained than the I = 
frequency shifts because of the small amplitude of the / = 3 
modes in the power spectrum. Besides, the I = 1 modes are 
observed to have a deeper downturn above 2100 //Hz than the 
I = modes. Nevertheless, these downturn and upturn, although 
consistent between the I = and I = 1 frequency shifts, occur 
when the SNR of the modes start to get small. Thus we cannot 
exclude any effects due to the low SNR to explain the behav- 
ior of the frequency shifts at high frequency. Longer datasets 



will be needed to improve the SNR of the measurements in or- 
der to confirm the observed shifts at high frequency. However, if 
they are real, these reversals of the sign at high frequencies can 
be qualitatively explained b y changes in chrom ospheric mag- 
netic field and temperature (iJohnston et aU[l995l) . The insert in 
Fig. [2] shows the increasing part of the mean frequency shift up 
to 2100 //Hz, which would reflect structural changes in and just 
below the photosphere with stellar a ctivity if we suppose si milar 
mechanisms as in the Sun (see, e.g jGoldreich et alTll991l) . The 
solid red line, which corresponds to a weighted linear fit, clearly 
shows the increase in frequency shift as a function of frequency. 
We note the excellent agreement with the results from the cross- 
correlation method (horizontal dashed lines). We also calculated 
the variations of the mean mode amplitude as a function of fre - 
quency using the same methodology as in iGarcia et al.l (|2010), 
but by separating the power spectrum into two regions. But no 
reliable results were obtained due to the very large uncertainties 
related to the reduced frequency ranges. 

Table Q] gives the mean frequency shifts obtained from the 
local analysis of the individual mode frequencies (method #1) 
averaged between 1470 //Hz and 1815 //Hz, and between 1815 
//Hz and 2160 //Hz, using the formal errors as weights. These 
frequency ranges were chosen to contain the same number of 
overtones of I = and I = 1 modes, and correspond to the 
represented frequency shifts in the insert of Fig. The re- 
sults obtained over the same frequency ranges from the analy- 
ses from the global fitting of the p-mode spectrum (method #2) 
and the cross-correlation technique (method #3) are also given 
in Table Q] The three independent methods agree within the 
lcr errorbars. The I = modes - and the even pair of modes 
(/ = 0,2) as well - are observed to have, in average, larger 
frequency shifts than the / = 1 modes. The total amount of 
change in mode frequency between periods of low and high stel- 
lar activity of HD 49933 is consistent with the temporal varia- 
tio n of the mode freq uencies with magnetic activity measured 
bv lGarcia et al.l (120101) . In order to check the consistency of the 
results, we also analyzed two subseries of 50 days each, centered 
on the period of low activity and on the rising phase of the mag- 
netic activity. Similar results as Fig.|2]and TableQ~]were obtained 
from the mode peak-fitting and the cross-correlation methods. 
For example, Fig. [3] compares the mean frequency shifts ob- 
tained from the analysis of the 90-day and 47-day subseries with 
those from the analysis of the two 50-day subseries in the case of 
the local mode peak-fitting analysis (method #1). The associated 
Pearson's linear correlation coefficient is about 0.89. 
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Fig. 2. Frequency shifts of the / = (blue squares) and / = 1 modes (red diamonds) of HD 49933 as a function of frequency 
obtained from method #1. The corresponding weighted mean values of the frequency shift between consecutive I = and 1 modes 
are represented by the black dots (solid line). A zoom-in on the frequency range up to 2100 //Hz is shown in the insert, which 
represents the mean frequency shift. The horizontal blue dashed lines correspond to the frequency shifts obtained with the cross- 
correlation method (#3). The solid red line corresponds to a weighted linear fit. 



Table 1. Mean frequency shifts (6vi) in //Hz over two frequency 
ranges obtained from the local and global mode peak-fitting, and 
the cross-correlation analyses for HD 49933. 



Method 


(Sv,) 


Frequency range 






OiHz) 






[1470-1815] 


[1815-2160] 


#1 


(I = 0, 1> 


0.26±0.35 


1.35+0.36 


#1 


</ = o> 


0.87±0.64 


2.26+0.61 


#1 


d=D 


0.01±0.42 


0.87±0.44 


#2 


(1 = 0,1,2) 


0.51±0.41 


1.23±0.43 


#2 


(even /> 


1.15+0.67 


2.34±0.82 


#2 


(odd /) 


0.12±0.53 


0.80+0.51 


#3 


(all I) 


0.19±0.12 


1.32+0.18 


#3 


(even /) 


0.64±0.17 


1.70±0.55 


#3 


(odd /> 


-0.05+0.15 


0.83±0.17 



Notes. The methods #1 and #2 correspond, respectively, to the local 
and global mode peak-fitting analyses (Sect. I2.lt . and the method #3 
corresponds to the cross-correlation analysis (Sect, l2~2~l . 
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Fig. 3. Comparison between the frequency shifts obtained from 
the analysis of the 90-day and 47-day subseries with those from 
the analysis of the two 50-day subseries in the case of method 
#1 . The solid line corresponds to the linear fit between the data, 
and the dashed line represents the 1 : 1 correlation. 



The frequency dependence of the frequency shifts measured 
in HD 49933 is comparable in shape to the one observed in 
the Sun, using both global (e.g.. Ilibbrecht & Woodardl l l990; 



Chao lin et alHl998tlGeUv et al.L l2002: Sala bert et al.Ll2004h and 



local (e.g., Ho we et al 1 120081) helioseismic measurements. That 
indicates then the presence of similar physical phenomena driv- 
ing the frequency shifts of the modes of oscillation as the ones 
taking place in the Sun, which are understood to arise from 
changes in the outer layers of the star along its activity cycle. 
However, the frequency shift measured in HD 49933 is at least 
five times larger than in the Sun, which is of about 0.5 //Hz at 
3700 //Hz between t he maximum and the minimum of the 11- 
year solar cycle (e.g. Sala bert et all 120041) . 



4. Summary 

The p-mode frequencies and amplitudes of the solar-like oscil- 
lating star HD 49933 were o bserved to vary in t ime with the 
magnetic activity of the star (IGarcia et all l2010h . We studied 
here the frequency dependence of th e p-mode frequency sh ift by 
analyzing the same observations as iBenomar et al 1 (12009b and 
IGarcia et al.1 (1201 Oh collected by the CoRoT satellite. The fre- 
quency shifts were obtained independently through the analysis 
of the mode frequencies extracted by both a local and a global 
peak-fitting procedures, and through a cross-correlation method. 
We showed that the frequency shifts in HD 49933 increase with 
frequency reaching a maximum of about 2 //Hz at 2100 //Hz. We 
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also inferred that the observed variations in mode frequencies 
are related to changes in the outer layers of HD 49933 due to its 
magnetic activity, as the activity-frequency relationship with fre- 
quency is analogous in shape to the one observed along the 11- 
year solar cycle. HD 49933 is thus the second star after the Sun 
for which the frequency dependence of the p-mode frequency 
shifts with magnetic activity has been measured. 
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